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It is shown that 21 cm intensity mapping could be used in the near term to make cosmologically 
useful measurements. Large scale structure could be detected using existing radio telescopes, or 
using prototypes for dedicated redshift survey telescopes. This would provide a measure of the 
mean neutral hydrogen density, using redshift space distortions to break the degeneracy with the 
linear bias. We find that with only 200 hours of observing time on the Green Bank Telescope, the 
neutral hydrogen density could be measured to 25% precision at redshift 0.54 < z < 1.09. This 
compares favourably to current measurements, uses independent techniques, and would settle the 
controversy over an important parameter which impacts galaxy formation studies. In addition, a 
4000 hour survey would allow for the detection of baryon acoustic oscillations, giving a cosmological 
distance measure at 3.5% precision. These observation time requirements could be greatly reduced 
with the construction of multiple pixel receivers. Similar results are possible using prototypes for 
dedicated cylindrical telescopes on month time scales, or SKA pathfinder aperture arrays on day 
time scales. Such measurements promise to improve our understanding of these quantities while 
beating a path for future generations of hydrogen surveys. 



I. INTRODUCTION 

An upcoming class of experiments propose the obser- 
vation of the 21cm spectral line over large volumes, de- 
tecting large scale structure in three dimensions [T] . This 
method is sensitive to a redshift range of z > 1, which 
is observationally difficult for optical experiments, due 
to a dearth of spectral lines in the atmospheric trans- 
parency window. Such experiments require only limited 
resolution to resolve the structures of primary cosmolog- 
ical interest, above the non- linear scale. At z w 1 this 
corresponds to tenths of degrees. There is no need to 
detect individual galaxies, and in general each pixel will 
contain many. This process is referred to as 21 cm inten- 
sity mapping. A first detection of large scale structure in 
the 21cm intensity field was reported in 2008 [2]. 

Intensity mapping is sensitive to the large scale power 
spectra in both the transverse and longitudinal direc- 
tions. From this signal, signatures such as the baryon 
acoustic oscillations (BAO), weak lensing and redshift 
space distortions (RSD) can be detected and used to gain 
cosmological insight. 

To perform such a survey dedicated cylindrical radio 
telescopes have been proposed [4], which could map 
a large fraction of the sky over a wide redshift range 
and on timescales of several years. These experiments 
are economical since their low resolution requirements 
imply limited size and they have no moving parts. It 
has been shown that BAO detections from 21cm inten- 
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sity mapping are powerful probes of dark energy, com- 
paring favourably with Dark Energy Task Force Stage 
IV projects within the figure of merit framework [5J [3] . 
Additionally, in [7] it was shown that such experiments 
could tightly constrain theories of modified gravity, mak- 
ing extensive use of weak lensing information. The f(R) 
theory, for instance, could be constrained nearly to the 
point where the chameleon mechanism masks any devia- 
tions from standard gravity before they first appear. 

While such experiments will be very powerful probes 
of the Universe, it is useful to explore how 21 cm inten- 
sity mapping could be employed in the short term. Here 
we discuss surveys that could be performed using existing 
radio telescopes, where the Green Bank Telescope will be 
used as an example. Prototypes for the above mentioned 
cylindrical telescopes (for which some infrastructure al- 
ready exists) are also considered. Finally, the Square 
Kilometre Array Design Studies (SKADS) focused on de- 
veloping aperture array technology for the Square Kilo- 
metre Array [5]. Motivated by the proposed Aperture 
Array Astronomical Imaging Verification (A3IV) project 
OHD], we consider pathfinders for such arrays. These 
aperture arrays could share many characteristics in com- 
mon with cylindrical telescope prototypes, but would be 
capable of much greater survey speed. 

While these limited resources would not have nearly 
the statistical power required to detect effects like weak 
lensing, a detection of the RSD would be possible. This 
would give a measure of the mean density of neutral hy- 
drogen in the Universe. This has been an important 
and controversial parameter in galaxy formation studies 
and a precise measurement would be invaluable in this 
field [11]. In addition BAO are considered, a detection 
of which would yield cosmologically useful information 
about the Universe's expansion history. 
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In this paper we first describe the RSD and BAO and 
the information that can be achieved with their detection. 
We then present forecasts for 21 cm redshift surveys as a 
function of telescope time, followed by a brief discussion 
of these results. 

We assume a fiducial ACDM cosmology with param- 
eters: Q m = 0.24, n b = 0.042, fl A = 0.76, h = 0.73, 
n s = 0.95 and log 10 A s — —8.64; where these represent 
the matter, baryon and vacuum energy densities (as frac- 
tions of the critical density), the dimensionless Hubble 
constant, spectral index and logarithm of the amplitude 
of scalar perturbations. 



II. REDSHIFT SPACE DISTORTIONS 

In spectroscopic surveys, radial distances are given by 
redshifts. However redshift does not map directly onto 
distance as matter also has peculiar velocities, which 
Doppler shifts incoming photons. On large scales, these 
velocities are coherent and result in additional apparent 
clustering of matter in redshift space. In linear theory, 
the net effect is an enhancement of power for Fourier 
modes with wave vectors along the line of sight [T2] , 



P* x (k,z) 



P(k,z) 



(1) 



Here P|- is the power spectrum of tracer X (assumed to 
be linearly biased) as observed in redshift space, P is the 
matter power spectrum and fit = k-z. The bias, 6, quan- 
tifies the degree to which the density perturbation of the 
tracer follows the density perturbation of the underlying 
dark matter. The redshift space distortion parameter /3 
is equal to f/b in linear theory, where / is the dimen- 
sionless linear growth rate. To a good approximation 

f( z ) = n m (z) - 55 . 

We define the signal neutral hydrogen power spectrum 

as 



P m (k, z) ee x 2 m py$,z), 



(2) 



where x^aiz) is the mean fraction of hydrogen that is neu- 
tral. It is the parameter xhi that we wish to measure. 
The above definition is useful because it defines the quan- 
tity that is most directly measured in a 21 cm redshift 
survey. Note that in discussions of the more standard 
galaxy redshift surveys one usually takes a measurement 
of the mean density for granted; however, 21cm inten- 
sity mapping will make differential measurements and 
will not measure the mean signal directly. One cannot 
divide out the mean to define a measurement of fluctua- 
tions around the mean. Using the fact that at late times, 
on large scales, structure undergoes scale independent 
growth, the above quantity can be written as 



= b 2 x 



g(z)a 
a 



l+p(4 P(k,z ), (3) 



where subscript refers to some early time well after 
recombination and g is the growth factor (relative to an 
Einstein-de Sitter Universe). 

Because the power spectrum at early times, P(k, Zq), 
can be inferred from the Cosmic Microwave Background 
(to good enough accuracy for our purposes in this pa- 
per), the observed power spectrum can be parameterized 
by just two redshift dependent numbers, /3 and the com- 
bination of scale independent prefactors in Equation [3j 

A H ee b 2 x 2 m g 2 . (4) 

It is these two parameters that can be determined from a 
21 cm redshift survey using the RSD. In general Ah will 
be better measured than j3 and does not significantly 
contribute to the uncertainty in xhi- 

The factors f{z) and g(z) depend on the expansion 
history. If one allows for a general expansion history (for 
instance, the WMAP allowed CDM model with arbitrary 
curvature and equation of state, OWCDM) these factors 
are poorly determined since there is currently little data 
that directly probes the expansion in this era. However, 
if one is willing to assume a flat ACDM expansion history, 
then the current uncertainty in the late time expansion is 
attenuated at the redshifts of interest (since dark energy 
is sub-dominant to matter at z = 1) and uncertainties in 
theses parameters can be ignored. 

As such, a measurement of /3 gives a measurement of 
the bias b which in turn gives a measurement of xhi ■ We 
have 

a (ACDM assumed). (5) 

xm P 

If one is not willing to assume an expansion history, it 
is a simple matter to propagate the corresponding uncer- 
tainties in the expansion parameters. 

To estimate errors on /3, we assume that the primor- 
dial power spectrum is essentially known from the cosmic 
microwave background, and we fix all parameters except 
for the amplitude A s , spectral index n s , and running of 
the spectral index a s . The observable power spectrum is 
then parameterized by /3, A s , n s and a s . We then use 
the Fisher matrix formalism to determine how precisely 
/3 can be measured from the 21 cm survey, marginalizing 
over the other three parameters and using no other in- 
formation. The parameter A s is used as a stand-in for 
the other parameters that affect the overall amplitude of 
the power spectrum: the bias and xhi- Its marginaliza- 
tion is critical to account for the fact that we have no a 
priori information about these parameters. The spectral 
index marginalization is not strictly necessary but allows 
for some degradation due to concern about scale depen- 
dence of bias. We have restricted the numerator of the 
Fisher matrix to include only the linear theory power as 
suppressed by the non-linear BAO erasure kernels of [T3] 
(all the linear power is included though, not BAO only), 
so non-linearity should not be a significant issue at the 
level of precision discussed here. This treatment of the 
non-linearity cutoff is also motivated by the propagator 
work of [HHS]. 
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FIG. 1: Baryon acoustic oscillations averaged over all direc- 
tions. To show the BAO we plot the ratio of the full mat- 
ter power spectrum to the wiggle- free power spectrum of [17] . 
The error bars represent projections of the sensitivity possible 
with 4000 hours observing time on GBT at 0.54 < z < 1.09. 



III. BARYON ACOUSTIC OSCILLATIONS 

Acoustic oscillations in the primordial photon-baryon 
plasma have ubiquitously left a distinctive imprint in the 
distribution of matter in the Universe today. This pro- 
cess is understood from first principles and gives a clean 
length scale in the Universe's large scale structure, largely 
free of systematic uncertainties, and calibrations. This 
can be used to measure the global cosmological expan- 
sion history through the angular diameter distance vs 
rcdshift relation. The detailed expansion will differ be- 
tween a pure cosmological constant and the various other 
cosmological models. 

We use essentially the method of [13] for estimating 
distance errors obtainable from a BAO measurement, 
including 50% reconstruction of nonlinear degradation 
of the BAO feature (although this is unimportant since 
experiments considered here have low resolution). The 
BAO feature is isolated by dividing the total power spec- 
trum |17j by the wiggle-free power spectrum and sub- 
tracting unity, as illustrated in Figure [I] The wiggles are 
then parameterized by an overall amplitude, and a length 
scale dilation (here A w and D respectively), which con- 
trol the vertical and horizontal stretch of the theoretical 
curve shown in Figure [I] Our errors on A w come from a 
straightforward extension of the [13 method for estimat- 
ing BAO errors. In addition to the BAO distance scale 
as a free parameter in our Fisher matrix, we include A w 
as a free parameter. This is similar to what one some- 
times tries to do by including the baryon/dark matter 
density ratio as a parameter, but more straightforward 
to interpret. 

The ability to measure A w (which is zero in the ab- 
sence of the BAO) represents the ability to detect the 
presence these wiggles. A measurement of D allows one 



to associate a comoving distance to length scales on the 
sky. This gives a measurement of the angular diameter 
distance {djCj for detections in the transverse direction, 
and the Hubble parameter (H) if the wiggles are detected 
in the longitudinal direction. 



IV. FORECASTS 

We present forecasts for the Green Bank Telescope and 
prototypes for two classes of telescope: cylindrical tele- 
scopes and SKA aperture arrays. The signal available for 
21 cm experiments is proportional to the neutral hydro- 
gen fraction and bias. For estimating telescope sensitivity 
we assume that the product of the bias and the neutral 
hydrogen density flmb = 0.0004 today[S], and that the 
neutral hydrogen fraction and bias do not evolve. These 
assumptions only affect the sensitivity of the telescopes 
and not the translation from uncertainty on Pgj to the 
uncertainty on xjji- Also, as in galaxy surveys, there is 
expected to stochastic shot noise component and we as- 
sume Poisson noise with an effective object number den- 
sity n = 0.03 per cubic /i _1 Mpc. Note that stochastic 
noise at this level is negligible, as should be the case in 
practice. 

The 21 cm intensity mapping technique is expected 
to be complicated by a variety of contaminating effects. 
These include diffuse foregrounds (predominantly galac- 
tic synchrotron), radio frequency interference and bright 
point sources. The degree to which these contaminants 
will limit future surveys has yet to be quantified, and 
here we simply ignore them. As such these forecasts are 
theoretical idealizations. Methods for dealing with these 
contaminants are discussed in [S]- 

The Green Bank Telescope is a 100 m diameter circu- 
lar telescope with a system temperature of 25 K. It has 
interchangeable single pixel receivers at the frequencies 
of interest with bandwidths of approximately 200 MHz. 
For extended surveys, multiple pixel receivers could be 
implemented. The construction of a four pixel receiver 
is within reason and would reduce the required telescope 
time by a factor of four. In planning a survey on GBT, 
it is important to choose an appropriate survey area. As 
illustrated in Figure [2j at fixed observing time, there is a 
survey area that best measures the desired parameters. 
For all results the survey area has been roughly optimized 
for the quantity being measured. The optimized areas are 
shown in Figure [3j Results are essentially insensitive to 
this area within a factor of 2 of the optimum. 

Prototyping for dedicated cylindrical telescopes is in 
its early stages. We present forecasts for a hypotheti- 
cal not-too-far-off telescope, composed of two cylinders. 
The total array measures 40mx40m with 300 dipole re- 
ceivers with 200 MHz bandwidth, and a system temper- 
ature of 100 K. Such telescopes have no moving parts 
and point solely by the rotation of the earth. As such, 
the area of the survey is set by latitude, receiver spac- 
ing, and obstruction by foregrounds; we assume 15 000 
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FIG. 2: Ability of GBT to measure the BAO and redshift 
space distortions as a function of survey area at fixed observ- 
ing time. Presented survey is between z = 0.54 and z = 1.09 
and observing time is 1440 hours. A factor of 10 has been 
removed from the A w curve. 
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FIG. 3: Roughly optimized survey area as a function of tele- 
scope time on GBT. Redshift range is between z — 0.54 and 
z = 1.09. 



square degrees. The survey area scales with the number 
of receivers leaving the noise per pixel unchanged. Thus 
the squared errors on measured quantities scale inversely 
with the number of receivers. Note that this is a slightly 
different scaling than the area optimized case of GBT, 
where it is the time axis that is scaled by the number of 
receivers. However, in the area optimized case variances 
also scale inversely with time as area optimization effec- 
tively fixes the noise per pixel. As such our results for 
GBT follow both scalings. 

The forecasts for cylindrical telescopes are also appli- 
cable to pathfinder aperture arrays [5]. Where-as cylin- 
ders use optics to form a beam in one dimension and in- 
terferometry in the other, aperture arrays directly sam- 
ple the incoming radio waves without optics. Interfer- 
ometry is used in both dimensions, forming a two di- 
mensional array of beams, instead of a two pixel wide 
strip for a two cylinder telescope. In principle it would 
be possible for an aperture array to monitor essentially 
the whole sky simultaneously, and thus form thousands 
of beams. In practise digitizing every antenna is costly 
and many antennas must be added in analogue in such 
a way that some beams are preserved but many are can- 
celled. We consider a compact aperture array that has 
the same area and resolution as the cylinder considered 
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FIG. 4: (Color online) Forecasts for fractional error on red- 
shift space distortion and baryon acoustic oscillation param- 
eters for intensity mapping surveys on the Green Bank Tele- 
scope (GBT). Frequency bins are approximately 200 MHz 
wide and correspond to available GBT receivers. Uncertain- 
ties on D should not be trusted unless the uncertainty on A w 
is less than 50% (see text). 



here. This is almost identical in scale as the proposed 
A3IV We assume preliminary A3IV specifications, 
with 700 effective receivers, 300 MHz bandwidth, and 
50 K system temperature [18] . The aperture array could 
thus perform the same survey as the cylinder but a fac- 
tor of (300 MHz/200 MHz)(700/300)(100K/50K) 2 = 14 
faster. Note however, that for aperture arrays there is 
added freedom in which beams are preserved when an- 
tennas are added. It would be thus possible to optimize 
the area of the survey as in the GBT case. 

Figures [4] and [5] show the obtainable fractional errors 
on RSD parameter f3 and BAO parameters A w , the am- 
plitude of wiggles, and D the overall dilation factor. D is 
defined as a simultaneous dilation in both the longitudi- 
nal and transverse directions, i.e. both H and cIa are pro- 
portional to D. There is another "skew" parameter which 
trades one for the other, such that the Hubble parameter 
and the angular diameter distance can be independently 
determined, however, this parameter is generally not as 
precisely measured [19 . D is the mode that contains 
most of the information available in the BAO and we 
marginalize over the skew parameter. The marginalized 
error on D is independent of the exact definition of the 
skew. Fractional errors on H and d,A are of order twice 
the fractional error on D. 

Referring to Figure [l] above, it can be seen intuitively 
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FIG. 5: (Color online) Forecasts for fractional error on red- 
shift space distortion and baryon acoustic oscillation parame- 
ters for intensity mapping surveys on a prototype cylindrical 
telescope. Frequency bins are 200 MHz wide corresponding to 
the capacity of the correlators which will likely be available. 
These result also apply to the aperture telescope but with 
the observing time reduced by a factor of 14. Uncertainties 
on D should not be trusted unless the uncertainty on A w is 
less than 50% (see text). Observing time does not account 
for lost time due to foreground obstruction. 

that the parameters D and A w are weakly correlated. 
Furthermore, since the signal is linear in the parameter 
A w , our linear Fisher analysis applies even for large er- 
rors. This however cannot be said about D. Indeed any 
uncertainty on D, that brings the smallest scale peak we 
resolve more than ~ n/2 out of phase, cannot be trusted. 
It can be seen from Figure [I] that this corresponds to a 
fractional error of order 10%, which of course depends 
on resolution. For this reason, we require that the un- 
certainty on A w must be at most 50% (a 95% confidence 
detection of the BAO) before we have any faith in the 
uncertainty in D. We note that Fisher analysis is not 
the optimal tool for determining when an effect can first 
be measured; however, it applies for any cosmologically 
useful measurements of the BAO. 



V. DISCUSSION 

Typically measurements of the neutral hydrogen den- 
sity at high redshift are made using damped Lyman-a 
(DLA) absorption lines; current measurement uncertain- 



ties being at the 25% statistical precision level in the 
0.5 < z < 2 redshift range (20] [5T]. However, it has been 
argued that these measurements are biased high [55] ren- 
dering the quantity effectively uncertain by a factor of 3. 
Hydrogen is the main baryonic component in the universe 
and it becomes neutral after falling into galaxies and be- 
coming self shielded from ionizing radiation. As such, the 
abundance of neutral hydrogen is linked to the availabil- 
ity of fuel for star formation [2"3l [5~4] . Understanding how 
the neutral hydrogen evolves over cosmic time is key to 
understanding the star formation history and feedback 
processes in galaxy formation studies [55] [5B] . Addition- 
ally, the linear bias, and any scale dependence it might 
have, gives valuable information about how the gas is 
distributed. Finally, these quantities are crucial for es- 
timating the sensitivity of future 21 cm redshift surveys 
since the signal is proportional to the product of the bias 
and the mean neutral density [5j. 

We have shown that even with existing telescopes, it 
is possible to use 21 cm intensity mapping to make useful 
measurements of large scale structure at high redshift. 
As seen in Figure [4] a 4rj detection of redshift space dis- 
tortions could be made at z ss 0.8 with only 200 hours of 
telescope time at GBT. This would provide a 25% mea- 
surement of the neutral hydrogen fraction in the Universe 
using methods independent of DLA absorption lines. A 
longer survey using 1000 hours of telescope time could 
make ^12% measurements. 

Surveys extending to this level of precision become cos- 
mologically useful. With 4000 observing hours on GBT 
(1000 hours with a four pixel receiver), the BAO overall 
distance scale could be measured to 3.5% precision at the 
same redshift. This is approaching the precision of the 
WiggleZ survey, which will make a ~2% measurement of 
this scale over a similar redshift range J27] [28 29 1 . Be- 
cause WiggleZ will make a similar measurement, such a 
survey would not have a dramatic effect on cosmologi- 
cal parameter estimations. However, it would provide an 
excellent verification of these measurements, using com- 
pletely different methods, in different regions of the sky, 
and at low cost. 

Prototypes for cylindrical telescopes could perform 
similar science to existing telescopes except — with dedi- 
cated resources — longer integration times would be fea- 
sible. This in part makes up for the limited resolution 
as for 40 m telescopes there is a substantial loss of in- 
formation, and only the first wiggles are resolved. The 
measurements described here would constitute a proving 
ground for this technique. The success of these proto- 
types would be a clear indicator of the power and future 
success of full scale cylindrical telescopes. 

The most powerful telescope considered is the aperture 
array, which would be capable of making sub present level 
BAO measurements with only a few weeks of dedicated 
observing. The design for demonstrator telescope A3IV 
has yet to be finalized but the proposed telescope is of 
nearly the same scale as the aperture array considered 
here. Depending on its eventual configuration, the A3IV 
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could be a powerful probe of the Universe. 

We have shown that 21 cm intensity mapping surveys 
could be employed in the short term to make useful mea- 
surements with large scale structure. With relatively 
small initial resource allocations, requisite techniques 
such as foreground subtraction can be tried and tested 
while performing valuable science. Such short term ap- 
plications of this promising method will lay the trail for 
future dark energy surveys. 



Acknowledgments 



We thank Ger de Bruyn for preliminary specifications 
of the A3IV. KM is supported by NSERC Canada. PM 
acknowledges support of the Beatrice D. Tremaine Fel- 
lowship. 



[1] J. B. Peterson et al. (2009), 0902.3091. [15 
[2] U.-L. Pen, L. Staveley-Smith, J. Peterson, and T.-C. 

Chang (2008), 0802.3239. [16 
[3] J. B. Peterson, K. Bandura, and U. L. Pen (2006), astro- 

ph/0606104. [17 
[4] H.-J. Seo et al. (2009), 0910.5007. 

[5] T.-C. Chang, U.-L. Pen, J. B. Peterson, and P. McDon- [18 
aid, Phys. Rev. Lett. 100, 091303 (2008), 0709.3672. [19 
[6] A. Albrecht et al. (2006), astro-ph/0609591. 
[7] K. W. Masui, F. Schmidt, U.-L. Pen, and P. McDonald, [20 

Phys. Rev. D81, 062001 (2010), 0911.3552. 
[8] A. Faulkner et al., Tech. Rep. DS8_T1, SKADS (2010), [21 

http : //www. skads-eu. org. 
[9] A. van Ardenne, in Wide Field Science and Technology [22 
for the SKA, edited by S. A. Torchinsky, A. van Ar- [23 
denne, T. van den Brink-Havinga, A. J. J. van Es, and 
A. J. Faulkner (Chateau de Limelette, Belgium, 2009), [24 
pp. 407-410. [25 
[10] Aperture array verification program, 

http://www.ska-aavp.eu/. [26 
[11] M. E. Putman et al. (2009), 0902.4717. [27 
[12] N. Kaiser, Mon. Not. R. ast. Soc. 227, 1 (1987). [28 
[13] H.-J. Seo and D. J. Eisenstein, Astrophys. J. 665, 14 [29 

(2007), arXiv:astro-ph/0701079. 
[14] M. Crocce and R. Scoccimarro, Phys. Rev. D 73, 063520 
(2006). 



M. Crocce and R. Scoccimarro, Phys. Rev. D 73, 063519 
(2006). 

M. Crocce and R. Scoccimarro, Phys. Rev. D 77, 023533 
(2008), arXiv:0704.2783. 

D. J. Eisenstein and W. Hu, Astrophys. J. 496, 605+ 
(1998). 

G. de Bruyn, private communication (2010). 

N. Padmanabhan and M. White, Phys. Rev. D 77, 

123540 (2008), 0804.0799. 

J. X. Prochaska, S. Herbert-Fort, and A. M. Wolfe, As- 
trophys. J. 635, 123 (2005), astro-ph/0508361. 
S. M. Rao, D. A. Turnshek, and D. B. Nestor, Astrophys. 
J. 636, 610 (2006), astro-ph/0509469. 
J. X. Prochaska and A. M. Wolfe (2008), 0811.2003. 
A. M. Wolfe, D. A. Turnshek, H. E. Smith, and R. D. 
Cohen, ApJS 61, 249 (1986). 

Y. C. Pei and S. M. Fall, Astrophys. J. 454, 69 (1995). 

A. M. Wolfe, E. Gawiser, and J. X. Prochaska, Ann. Rev. 

Astron. Astrophys. 43, 861 (2005), astro-ph/0509481. 

S. Shen, J. Wadsley, and G. Stinson (2009), 0910.5956. 

C. Blake et al. (2009), 0901.2587. 

M. J. Drinkwater et al. (2009), 0911.4246. 

The projected uncertainties that include reconstruction 

in these references are on H and dA ■ The uncertainty on 

D is inferred from these. 



